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ABSTRACT 

Context. If one wants to understand the physics of blazars, better simultaneous observations are important at all wavelengths, so it 
was fortunate that a ToO observation of the TeV-emitting blazar Mrk 421 with INTEGRAL could be triggered in June 2006 by an 
increase in the RXTE count rate to more than 30 mCrab. The source was then observed with all INTEGRAL instruments, with the 
exception of the spectrometer SPI, for a total exposure of 829 ks. During this time several outbursts were observed by IBIS and JEM- 
X. Multiwavelength observations were immediately triggered, and the source was observed at radio, optical, and X-ray wavelengths 
up to TeV energies. 

Aims. The data obtained during these observations were analysed with respect to time variability, time lags, correlated variability, and 
spectral evolution and then compiled in a v F^spectrum. 

Methods. The observations of the different instruments/telescopes were analysed with the usual correlation and time-analysis methods. 
The spectral analysis of the X-ray data was performed with XSPEC. 

Results. Four strong flares at X-rays were observed that were not seen at other wavelengths (partially because of missing data). 
From the fastest rise in the X-rays, an upper limit could be derived to the extension of the emission region. A time lag between 
high-energy and low-energy X-rays was observed, which allowed an estimation of the magnetic-field strength. The spectral analysis 
of the X-rays revealed a slight spectral hardening of the low-energy (3 — 43 keV) spectral index. The hardness-ratio analysis of the 
Swift-XRT (0.2 - 10 keV) data indicated a small correlation with the intensity; i. e., a hard-to-soft evolution was observed. At the 
energies of IBIS/ISGRI (20 - 150 keV), such correlations are less obvious. A multiwavelength spectrum was composed and the X-ray 
and bolometric luminosities calculated. 

Conclusions. The observed flaring activity of Mrk 421 is mainly visible at X-rays. It is found that the spectral change with intensity 
is small. But at least one flare showed a completely different spectral behaviour than the other flares, so one can conclude that each 
blob of relativistic-moving electrons has its own individual physical environment that leads to different emission characteristics. From 
a fit of a leptonic emission model to the data, one finds that the observed variability may be due to a varying efficiency of particle 
acceleration. 
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1. Introduction 

In 1992 the first AGN, Mrk 421, was detected at energies 
>500 GeV (Punch et al. [1992]). Since then, 20 detections of 
AGNs have been reported at TeV energies (Blazejowski et al. 
120051 Padovani et al. 120571 1 Wagner 200"7) . Seventeen of these 
objects belong to the HBL Lac type of AGNs. They are radio- 
loud sources with the radio emission originating mainly in a core 
region and they are characterised by a high polarization at ra- 
dio and optical wavelengths and a strong variability at all wave- 
lengths. The spectral characteristics point to non-thermal emis- 
sion processes that presumably take place in a narrow relativis- 
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tic jet pointed at a small angle to the line of sight. The spectral 
energy-density distribution of these sources shows two smooth 
broadband emission components: a first one that reaches a broad 
peak in the IR to X-ray region and a second one at GeV to TeV 
energies. Both emission components are clearly separated (see 
Figure 3 of lGhiselhni et al. 19961 ). 

It is believed by many researchers that both components 
are generated by the same leptonic population (electrons and/or 
positrons) moving at relativistic speed in the jet (Blandford and 
Rees 1978), which creates the low-energy photons via incoher- 
ent synchrotron radiation and the high-energy gamma-rays via 
inverse-Compton scattering of soft photons. This common origin 
would explain the similarity of the two components (Ghisellini 
and Maraschi 1996); hence, a similar temporal evolution is ex- 
pected for both components if this scenario is correct. Such 
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time correlations were indeed observed by observations at X- 
ray and TeV energies (Buckley et al.[l996; Catanese et al. 1 19971 
Foss ati et al. 20 08 ), and they even extend to the optical and ra- 
dio wavelengths (Katarzynski et al. 120031 1. This idea is further 
supported by the observed polarization at radio wavelengths 
dPiner and Edw ards 2005 ). However, there are some exceptions 
to these clear correlations, such as the orphan TeV flares reported 
by Blazejowski et al. (2005J and |Krawczynski et aL~T2 004) 
Unfortunately, measurements between 10 keV and 1 GeV that 
would allow us to further constrain this common origin have 
been scarce. Until now only few measurements between 10 keV 
and 1 GeV exist. COMPTEL has detected Mrk 421 in the 10-30 
MeV range with 3.2cr (Collmar et al. [19991 . EGRET, however, 
has detected this source many times above 100 MeV (Hartman 
et al. [19991 . 

Since the high-energy emission is not as well understood, 
many different models try to explain this emission component. 
In one scenario the low-energy synchrotron photons are boosted 
to high energies by the same electron population that creates 
the synchrotron photons [the synchrotron self-Compton (SSC) 
mode ls (Konigl [19811 IMaraschi et all 1 999, Bloom & Marscher 
1996)]. Although this process must be at work in all blazars 
and if the synchrotron-emission hypothesis is correct, it may 
not be the dominant one. In another model the seed photons for 
the inverse-Compton effect enter from outside the jet region, e. 
g. from the accretion disc or from clouds surrounding the jet 
(Dermer et al. 1992; Sikora et al. 1994). Apart from these lepton 
models, the so-called hadron models were proposed in which 
the high-energy y-rays are produced by proton-initiated cas- 
cades (Mannheim l 1 9931 Miicke & Protheroe 2000]( and/or proto- 
synchrotron emission (lAharonian 2 000 ; Mii cke et al. 2 003 ). 

The homogeneous SSC model makes very definite predic- 
tions about the correlated behaviour of the high-energy end of 
both the synchrotron and SSC components: simultaneous vari- 
ability of photons and well-defined correlated spectral changes 
in the medium/hard X-rays and TeV bands. But these models 
have problems with the bulk-Lorentz factor statistics (Henri & 
Sauge 2006) and it seems to be very important to measure a de- 
tailed variability pattern from X-rays up to TeV y-rays. Several 
alternatives have been proposed, and each needs to be quanti- 
fied in more detail. However, in general one can say that in the 
leptonic models synchrotron radiation is "primary" and inverse 
Compton radiation "secondary", while for at least one hadronic 
model it is the opposite: the y-rays are produced by the "pri- 
mary" protons through cascades, while the X-rays have a syn- 
chrotron origin in the "primary" electrons. In this case the rela- 
tion between the two components can be looser than for lepton 
models. Also, in this case the "cascade" spectra are quite soft 
(lots of soft y-rays). 

To investigate all these theoretical ideas and models, an 
INTEGRAL ToO proposal was submitted to observe the blazar 
Mrk 421 when it becomes active. Mrk 421 has a redshift of 
0.031, and with a distance of about 125 Mpc [for Ho = 71 
km/(s Mpc)] IS one of the closest and therefore brightest 
blazars. It hosts a supermassive black hole with a mass of (2-8) 
•10 8 M Q in its c entre dBarth et al. 2003] IFalomo et al. 20021 
ITreves et al. 20031) . At high energies, MRK 421 has been 
detected in the 10-30 MeV range by COMPTEL with 3.2cr 
(Collmar et al. 119991 ) and many times above 100 MeV by 
EGRET (Hartman et al [19991 and in th e TeV range by 
Cherenkov telescopes ( Kerrick et al. 19951 lAharonian et al.l 
1999; ISambruna et al. 2000[ IKrennrich et al. 200 11 
lAharonian et al. 20051 and Rebillot et al. 120061 ). In June 
2006, Mrk 421 became active and our INTEGRAL proposal 



was activated. The proposal foresaw not only observations with 
the INTEGRAL instruments but also observations with other 
telescopes at all wavelengths. Data at radio, optical, X-ray, 
and TeV wavelengths were taken.However, the observations 
at TeV energies turned out to be sparse because the visibility 
period for the Cherenkov telescopes approached its end. It was 
thus impossible to investigate possible correlations or time lags 
between the synchrotron and inverse-Compton emission as we 
wished. We therefore had to concentrate our analysis mainly 
on the INTEGRAL data itself and on the multi-wavelength 
spectrum. Preliminary results of this analysis have been already 
published by |Lichti et al. (2007)| 

2. OBSERVATIONS 

In April 2006 the blazar Mrk 421 increased its intensity to a 
level >30 mCrab as measured by the all-sky monitor (ASM) of 
RXTE (the quiescent intensity fluctuates strongly around an av- 
erage value of 10 - 15 mCrab). It remained at this level until 
September 2006. This triggered an INTEGRAL observation and 
correlated multiwavelength observations in the radio, optical, X- 
ray and TeV ranges. The details of the various observations are 
given in this section. 

2.1. INTEGRAL observations 

On June 14, 2006 (MJD 53900) at 09:53:16.8 hours UT, an ob- 
servation of this source was triggered with INTEGRAL. The 
observation lasted about nine days and ended on June 26, 2006 
(MJD 5391 1) at 02:03:55.4 hours UT. In total the source was ob- 
served for 829 ks with the instruments IBIS, JEM-X, and OMC 
in pointing mode (the spectrometer SPI was disabled because of 
annealing its Ge detectors). 

The data of the three operating instruments were analysed 
with the INTEGRAL Off-line Scientific Analysis (OSA) soft- 
ware version 5.1 using the latest response matrices available for 
that software. The collected data were screened by computing 
the median count rate for each science window, and it was then 
compared with each of the other science windows and their dis- 
tributions. Those science windows showing rates higher than 10 
standard deviations from the median count rate were checked 
again. By this screening effort, three suspicious science windows 
were detected. A closer look revealed that these science windows 
belonged to times when INTEGRAL was close to the radiation 
belts and that the high counting rates were due to the trapped 
charged particles. After this screening process, 230 science win- 
dows remained and were considered for further analysis. 

Mrk 421 was clearly detected by all three INTEGRAL in- 
struments. In Figure[T]IBIS/ISGRI, JEM-X, and OMC skymaps 
in galactic coordinates are shown. In all three maps the source 
is detected at a high significance level [by IBIS in the energy in- 
terval 20-40 keV with 160cr, by JEM-X in the energy interval 
3-4.5 keV with 720cr, and by OMC (limiting magnitude «18) 
with m v * 12.87)]. It should be mentioned that in the field of 
view of IBIS another AGN, NGC 4151, was observeded, which 
was detected in the energy intervals 20-50 keV, 50-100 keV, and 
100-150 keV with 34.7<x, 15.1cr, and 5.4<x, respectively. 

The magnitudes m of OMC were converted into differen- 
tial flux values with the following formula (obtained from A. 
D. Garau of the OMC team) 



S[ ^r] = 3.64 ■ 10~ 9 • 10- a4m (1) 

cm 2 sA 
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AS = 0.921 S -Am. (2) 

The calibration of (1) was obtained from the flux of a star with 
m = at 5500 A from lWamsteker 19811 The contribution of the 
host galaxy was estimated from the colour of a typical elliptical 
galaxy (V-R = 0.6) and its R-band flux of 16.5 mJy to be 9.5 
mJy. This host-galaxy flux was taken into account when calcu- 
lating the energy-density flux. A correction for a possible contri- 
bution from the bright (V = 6) close-by (~ 2 arcmin away) star 
5 1 Uma (HR 4309) is not necessary because of the sharp point- 
spread function of the OMC, which drops practically to zero at 
a distance of 2 arcmin (attenuation factor < 10~ 23 ). 

Immediately after the beginning of the observation with 
INTEGRAL, the other observers participating in this multi wave- 
length campaign were alerted, and they commenced with their 
observations as allowed by the observational constraints. 

2.2. VLBA observations 

Following the notification of the INTEGRAL observations, a re- 
quest to trigger our previously-approved target of opportunity 
VLBA observations was sent to NRAO on June 13, 2006. After 
preparation of the observing schedule, the project entered the 
queue for dynamical scheduling at the VLBA and actually got 
scheduled on June 25. The observations started at 17:00 UT on 
that day and lasted for 14 hours, with scans alternated on six dif- 
ferent frequency bands (2.3, 5.0, 8.4, 15.4, 22.2, and 43.2 GHz). 
The total observing time on each frequency band was roughly 
1 hour at 2.3, 5 and 8.4 GHz, 2 hours at 15 GHz, 3 hours at 
22 GHz, and 6 hours at 43 GHz. Additional scans on calibrator 
sources were also inserted into the schedule for proper amplitude 
and polarization calibration. 

After correlation at the Array Operation Center in Socorro, 
the data were calibrated and corrected for residual delay and 
delay rate using the NRAO Astronomical Image Processing 
System (AIPS) and subsequently imaged using the Caltech 
VLBI imaging software DIFMAP. The flux density was then 
integrated over each image to produce the total flux densities 
used in the present study of the spectral energy distribution of 
Mrk 42 1 . The values of the flux density derived so far in this way 
are 0.283 Jy at 5 GHz and 0.273 Jy at 8.4 GHz. Presentation of 
the VLBI images is deferred to a future paper, which will present 
and discuss the results of the VLBI data analysis in more details. 

2.3. Observations with the Metsahovi radiotelescope 

The 37 GHz observations were made with the 13.7 m diam- 
eter Metsahovi radio telescope, which is a radome-enclosed 
paraboloid antenna situated in Finland (24 23' 38"E, +60 13' 
05"). The measurements were made with a 1 GHz-band dual- 
beam receiver centred at 36.8 GHz. The HEMPT (high-electron 
mobility pseudomorphic transistor) front end operates at room 
temperature. The observations are ON-ON observations, alter- 
nating the source and the sky in each feed horn. A typical inte- 
gration time to obtain one flux-density data point is 1200-1400 s. 
The detection limit of our telescope at 37 GHz is approximately 
0.2 Jy under optimal conditions. Data points with a signal-to- 
noise ratio < 4 were handled as non-detections. 

The flux-density scale was set by observations of the bright 
compact HII region DR 21. Sources 3C 84 and 3C 274 were 
used as secondary calibrators. A detailed description on the data 
reduction and analysis is given in Terasranta et al. |(1998)[ The 
error estimate in the flux density includes the contribution from 



the measurement rms and the uncertainty of the absolute calibra- 
tion. 

The observations of Mrk 421 during the INTEGRAL cam- 
paign were made in varying weather conditions. The first mea- 
surements were non-detections, but shortly the flux level in- 
creased and several detections were achieved on days 53902 
and 53903. The rest of the campaign was lost due to changeable 
weather and the end of the telescope time allocation. 

The sparse radio data do not allow any conclusion by them- 
selves. For the Metsahovi observations Mark 421 stays basi- 
cally on the verge of the detection limit and it never gets very 
bright. This is typical of the TeV blazars observed at Metsahovi. 
However, this kind of behaviour is different from other GeV- 
peaking gamma-ray blazars, and may reflect the different mech- 
anism for the production of gamma rays in Mrk 42 1 and the other 
TeV blazars. 

2.4. KVA observations 

Mrk 421 was observed in the R-band (centred at 640 nm) with 
the 35 cm KVA telescope on La Palma, Canary Islands, dur- 
ing 8 nights on June 14-22, 2006. The camera employed was 
an SBIG ST-8 with a gain factor of 2.3 e /ADU and and read- 
out noise of 14 e~. Each night 4-10 exposures of 180 s were 
made. The images were processed in a standard way (bias sub- 
traction, dark subtraction and flat-fielding) and the raw counts 
of Mrk 421 and stars 1-3 in Villata et al. |(1998)| were integrated 
within a 15 arcsec diameter aperture. The magnitude of Mrk 421 
was determined by comparing its counts to star 1 , for which R = 
14.04 from Villata et al. (1998) was assumed. The use of differ- 
ential mode effectively eliminates the influence of varying trans- 
parency due to clouds, etc., and accurate fluxes can even be ob- 
tained under varying observing conditions. The flux of Mrk 421 
was then computed from 

F[Jy] = 3080.0 • 10" a4s , (3) 

where R is the R-band magnitude (with an effective wavelength 
of 6400 A). The normalization constant was taken from Table 
IV of Bessell (1979) The host galaxy flux with the 15 arcsec 
diameter aperture is 8.0+0.5 mJy (Nilsson et al. 2007), and this 
flux was subtracted from the observed fluxes before proceeding 
with the analysis. 

2.5. RXTE observations 

The data of the all-sky monitor (ASM) of RXTE (for a descrip- 
tion see Levine et al. 1 19961 1 monitors the X-ray sky more or less 
continuously in the energy range 1.5-12 keV. Its data are pub- 
licly available. The data for Mrk 421 were downloaded and the 
fifth lightcurve at the top of Figure [2] was produced. The RXTE 
counting rates were translated into energy fluxes by using the 
Crab-nebula flux from Kirsch et al. |(2005)[ These authors de- 
rived an energy flux of 2.2 ■ 10 -8 erg/(cm z s) for the Crab for the 
energy interval 2-10 keV (i. e. the energy range of the ASM). 
This corresponds to an RXTE counting rate of 73 counts/s, so 
1 RXTE count/s corresponds to 3 • 10~ ll) erg/(cm 2 s). With this 
value, the RXTE-counting rates were converted to energy fluxes. 

2.6. Swift-XRT observations 

The X-Ray Telescope (XRT) instrument on the Swift 
Observatory (Burr ows et al. 20031 and IGehrels et al. 20041 is 
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Fig. 1. IBIS/ISGRI, JEM-X, and OMC skymaps (from left to right) of the region around Mrk 421 for 20-40 and 3-4.5 keV and 
at visible wavelengths, respectively, in galactic coordinates. In the OMC map, the bright star 51 Uma (V = 6), which is clearly 
separated from Mrk 421 (distance ~2 arcmin), is seen at the edge of the field of view. 



sensitive to X-rays in the 0.2-10 keV band. These data were re- 
duced using the latest HEAsoft tools (version 6.1.0), including 
Swift software version 2.0 and the latest response (version 8) 
and ancillary response files (created using xrtmkarf) available 
in CALDB at the time of analysis. Data were screened with 
standard parameters, including the elimination of time periods 
when the CCD temperature was warmer than -48° C. Only WT 
(windowed timing) mode data was used in this analysis due to 
the high rates of the active source, and only grades 0-2 were 
included. Since the count rate stayed below ^100 c/s, the WT 
mode data is free of significant pile-up effects. These data were 
corrected for effects due to bad columns and bad pixels. Source 
and background regions were both chosen in a way that avoids 
overlap with serendipitous sources in the image. All analysis and 
fitting of XRT data was done in the 0.3 to 10 keV energy band. 
Due to the low hydrogen-column density towards Mrk 421 of 
1.43 • 10 20 cm~ 2 , the absorption is <0.5% and was neglected. 
Finally the measured counting rates were transformed into en- 
ergy fluxes using the formula 



F\ 



erg 
(cm 2 s) 



] = (4 + 2.5) • 1(T 



■R. 



cts/s- 



(4) 



This conversion factor was calculated using a subset of these 
Swift-XRT Mrk 421 observations, which included both highl- 
and low-state flux measurements. For each time period used, the 
rate was calculated using the analysis described above, and the 
spectral fits to these data allowed the flux to be calculated, thus 
resulting in a conversion factor. Since a single conversion fac- 
tor was used, rather than many individual time-resolved spectral 
fits, the large error bars are required to represent the variation in 
this conversion factor due to spectral variability of the source. 

2.7. Swift-BAT observations 



The Burst-Alert Telescope (BAT |Barthelmy et al. 2 005) 



on board the Swift satellite mission ( Geh rels et al. 20041 > is a 
coded-mask telescope sensitive to the 15-200 keV energy range. 
Thanks to the pointing strategy and to its wide field of view 
(FoV), BAT surveys ~80% of the sky every day. We thus looked 
in the archive for BAT observations that contained Mrk 421 in 
the FoV. We selected all observations included in the time-span 
June 12-26, 2006. The BAT data were processed using the 
HEASOFT 6.2 package and according to the recipes presented in 
http://swift.gsfc.nasa.gov/docs/swift/analysis/threads/bat_thread 
s.html. Spectra and lightcurves were corrected for off-axis 



Table 1. The measured TeV intensities (in Crab units) of Mrk 
421 as measured by the Whipple observatory. 



MJD^ 


flux value 


error of 


energy density 


energy-den- 




(Crab units) 


flux value 


[erg/(cm 2 s)] 


sity error 


53904.17 


0.25 


0.38 


2.90- 10-" 


4.46 • 10-" 


53905.17 


1.63 


0.42 


1.89- lO" 10 


6.31 • 10-" 


53907.20 


0.39 


0.28 


3.71 • 10-" 


2.82 ■ 10-" 



" The integration time for each data point was 28 minutes. The centre 
of each time interval is given. 



variation of the rates and residual background contamination as 



described in Ajello et al. (2007) 



2.8. Whipple observations 

The Whipple Telescope is located at the Fred Lawrence Whipple 
Observatory, at an elevation of 2300 m in Southern Arizona. 
It comprises a 10 m dish on which 248 mirrors are mounted. 
These mirrors reflect the Cherenkov light from extensive air 
showers onto a 379-element imaging camera at the focal plane 
of the telescope. The instrument is described in detail by 
IKildeaetal. (2007)1 

The Mrk 421 data presented here were taken in "tracking" 
mode at large zenith angles on June 18, 19, and 21, 2006 (MJD 
53904, 53905, and 53907). The observing scans were performed 
with the gamma-ray source at the centre of the field of view. 
No separate control data were taken. Rather, the background 
rate of gamma-ray-like events was estimated from the distri- 
bution of events passing all but the orientation image-selection 
cuts |(Horan et al. 20 02) The gamma-ray rates are calculated in 
units of the flux from large zenith-angle observations of the Crab 
Nebula, the standard candle in TeV astronomy, and are given in 
Table [1] They were converted into energy-flux values by inte- 
grating the Crab spectrum of Hilla s et al. (199 8) from the thresh- 
old energies of 0.6 and 0.9 TeV to infinity and using the result- 
ing flux values as normalization for the time-averaged spectrum 
of Mrk 421 as given by Ah aronian et al.l (1999). Under the as- 
sumption that the spectral shape is constant this normalization 
was adapted to the measured flux values. 

3. Results of the timing analysis 

It is known that blazars are time variable on all different time 
scales. This is especially true for TeV blazars and thus also for 
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Table 2. Weighted means of the fluxes/magnitude/counting rates 
in Figure [2] 



instrument 


average value 


error 


radio 


0.33 Jy 
1.83 • 10~ 2 Jy 


0.026 Jy 


KVA 


3 • 10- 4 Jy 


OMC 


2.579-10- 14 erg/(cm 2 sA) 


0.007 erg/(cm 2 s A) 


RXTE 


2.86 cts/s 


0.3 cts/s 



Mrk 421. We have therefore investigated our data in this respect 
and compare the data collected at different energies with each 
other and search for possible correlations betweeen them. The 
results of this exercise are presented in this section. 

3.1. The Lightcurves 

The lightcurves of the various observations are shown in 
Figures [2] to [7] Whereas the sampling of the lightcurves of the 
three INTEGRAL instruments and the SWIFT-XRT telescope is 
ample, this is not the case for the other lightcurves. This makes 
the search for intensity correlations between the different wave- 
bands difficult if not even impossible. Another fact is striking 
when one looks at the lightcurves. Although four strong flares 
are seen in the ISGRI lightcurve (and even more in the Swift- 
XRT lightcurve of Figure [3}, these flares are not visible at the 
other energies. An inspection of the X-ray lightcurves reveals 
the following: 

- the strength and the shapes of the flares vary significantly; 

- the shape of the flares is not symmetric (some flares show a 
sharp rise and a slow decay time, others a slow rise and a fast 
decay time); 

- the duration of the flares is about 0.5 days; 

- comparing Figures[3]and[4]one can see that the energy spec- 
trum of the flare around day 53910.5 is qualitatively softer 
than the spectrum of the other flares; 

- from a linear fit to all the data, it becomes obvious that the 
total intensity is slightly increasing with time. 

Especially interesting is that these flares are also not seen 
by the OMC. It is the common understanding that in blazars, if 
viewed under small angles, the optical thermal emission of the 
accretion disk and of the galaxy is outshined by the non-thermal 
beamed emission of the jet (Tavecchio 120051 1. Since both the 
X-rays and the optical photons are produced by the same pop- 
ulation of relativistic electrons via synchrotron emission, one 
would expect similar (if not even identical) lightcurves, if this 
scenario is correct. For Mrk 421, this happened to be the case 
in 2001 when a flare showing variability in the visible band co- 
incident with that at radio wavelengths and at TeV energies was 
observed (Katarzynski et al. 2003). That this is not observed in 
our data contradicts to this hypothesis and asks for an explana- 
tion. Possible explanations will be discussed in the last section. 

From the radio, KVA, OMC, and RXTE lightcurves of 
Figure [2] the average fluxes/magnitude/counting rates were cal- 
culated. The weighted means are given in Table [2] The ISGRI 
lightcurve of Figure|5]was used to define by eye the times when 
the source was in quiescent and active states. Three time inter- 
vals were specified and their boundaries are given in Table [3] 

3.2. Rise-time analysis 

The rise-time scale gives information about the extension of 
the emission region. From observations at TeV energies it is 
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Fig. 3. Swift-XRT lightcurve in the energy interval 0.3 - 10 keV. 
The dotted lines give only the coordinates, not the intervals of 
the active phases as in Figures [4]to[6] 
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Fig. 4. JEM-X lightcurve in the energy interval 10-15 keV. The 
time intervals for the active phases are indicated by the dotted 
lines. A parabola was fitted to all counting rates to guide the 
eyes. 
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Fig. 5. IBIS/ISGRI lightcurve in the energy interval 20-40 keV. 
The time intervals for the active phases are indicated by the dot- 
ted lines. A parabola was fitted to all counting-rate values to 
guide the eyes. 
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Fig. 2. All lightcurves for Mrk 421 at radio, optical, X-ray, and gamma-ray wavelengths. The data of the OMC [given in erg/(cm 2 s 
A~)l and KVA telescones were not extinction-corrected, since the extinction E»_v in the direction of Mrk 421 is onlv 0.03. The time 
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Fig. 6. IBIS/ISGRI lightcurve in the energy interval 40-100 keV. 
The time intervals for the active phases are indicated by the dot- 
ted lines. A parabola was fitted to all counting-rate values to 
guide the eyes. 
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Fig. 7. Swift-BAT lightcurve in the energy interval 15 - 200 keV. 

Table 3. Time boundaries in modified Julian Date (days:hours: 
minutes:seconds) for the active state. 



interval number 


start time 


end time 




1 


53902:07:47:23 


53904:00:29 


59 


2 


53906:09:53:24 


53906:19:07 


15 


3 


53909:09:36:14 


53910:18:49 


56 



known that Mrk 421 shows short-time variability with time 
scales of about one day and shorter (Aharonian et al. 2002, 
120031 ). Doubling times with time scales of about one day or 
even down to 15-20 minutes were also reported by Cortina and 
Schweizer p002)l and Krennrich et al. |(2003)| In May 1996 two 
fast flares were observed (Gaidos et al. 1996). The first one had 
an intensity 10 times higher than the one of the Crab nebula (the 
highest flux ever observed from this source at these energies) 
with a time scale of < 1 day. The second one a week later had an 
unprecedented time scale of ~30 minutes indicating a very small 
emission region! 

An inspection of the ISGRI ligtcurves reveals that the steep- 
est ascent occurred in the third time interval. Fitting an expo- 
nential rise-time law of the form a - e 1 ^ to the data of the energy 
interval 40 - 100 keV one finds a value of 2.79 hours for to. In 
the energy interval 20 - 40 keV, the increase is somewhat slower 



Table 4. The fractional variability measure V for the different 
lightcurves. 



instrument 


low-energy 
bound [eV] 


high-energy 
bound [eV] 


V 


Metsahovi 






undefined 


KVA 


1.85 


2.04 


0.0314 


OMC 


2.14 


2.39 


0.0393 


Swift-XRT 


300 


10000 


0.266 


RXTE 


2000 


10000 


0.4258 


JEM-X 


10000 


15000 


0.4266 


IBIS/ISGRI 


20000 


40000 


0.5146 


IBIS/ISGRI 


40000 


100000 


0.5723 



(only 4.28 hours). One can now calculate the size of the emission 
region with the formula 



/ < 



c ■ t o • 6 
1 +z 



(5) 



with c the speed of light, 6 the Doppler factor (=12; from Kino 
et al. 120021 ) and z the redshift. Inserting the value from above 
one means that the emission region must be < 234 AU (the 
Schwarzschild radius of Mrk 421 is between 4 and 16 AU and 
the radius of the last stable orbit between 12 and 47 AU). From 
the shortest time scale measured at TeV energies an extension 
was derived of the emission region of < 2.5 AU (Aharonian et 
al. 2002). From pure high-resolution radio-interferometric ob- 
servations Chariot et al. (2006) derived an upper limit on the 
high-energy emission region of 0.1 pc (=; 20000 AU). Thus a de- 
pendence of the size of the emission region on energy seems to 
exist: the higher the energy, the smaller the emission region! 

3.3. Variability analysis 



The fractional variability V (see appendix in Fossati et al. 2000 
for the definition) was calculated for all lightcurves. The result 
is shown in Table [4] The same behaviour as already found by 
previous authors (e. g. Fossati et al. 12000 Sembay et al. 2002) 
is confirmed by our observations: the fractional variability in- 
creases with energy. This behaviour seems to be common in 
blazars (Ulrich et al. 119971 ). If we fit a power law to the data 
we get V oc E 28 . This is in good agreement with the laws found 
by Fossati et al. J2000b (V oc E 025 ) and by Giebels et al. J2007b 
(V oc E 024 ). This reflects that the higher energies are emitted 
from regions closer to the central black hole, which are smaller 
in extension than regions farther away. Since the emission from 
small regions can vary faster than the ones from larger regions 
the observed trend can be understood qualitatively. 

3.4. Time-lag analysis 

The time lag between the two ISGRI lightcurves of the energy 
intervals 20-40 keV and 40-100 keV was calculated with the 
Z-transformed discrete-correlation function (ZDCF; Edelson & 
KrolikQim Tal l2"007i ). The value of the ZDCF as a function of 
the time lag t is shown in Figure [8] Five maxima with ZDF val- 
ues around 0.5 are seen, two at r « + 3.4 days, two at r ~ + 
7 days, and one at t « day. The shifts at + 3.4 days and 
± 7 days correspond to the time difference between the flares 
seen in Figures [5] and [6] and thus do not reveal new information. 
Important, however, is the question whether the central peak has 
its maximum really at t ~ days. To find this out, a Gaussian 
and a linear background were fitted to the data from r = — 1 day 
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Fig. 8. The Z-discrete correlation function of the IBIS data. 



to t = + 1 day. The data can be described well with a Gaussian 
ix 2 red = 0.114). The fit to the data reveals that, with a signifi- 
cance of 1 .8cr, the maximum is not at a time lag of zero, but at 
r = (-40.8 + 23.3) minutes. The negative sign indicates that the 
variations at lower energies lag behind those at the high energies. 

The same time-lag analysis was performed for the IBIS (20- 
40 keV) / JEM-X (10-15 keV) and the JEM-X (10-15 keV) / 
Swift-XRT (0.3-10 keV) lightcurves. In the first case the time 
lag is consistent with a lag of zero minutes (actually -4.25 + 
30.4 minutes), and in the second case the time lag is positive 
and has a value of (102.6 + 59) minutes. Similar time lags 
(35 - 47 minutes) were found by Fossati et al. (2000) for the 
energy ranges 0.1-1.5 keV and 3.5-10 keV using BeppoSAX 
data and by Takah ashi et al. (1996)| as a function of energy for 
the energy range 0.5-7.5 keV using ASCA data. Especially 
Takahashi et al. (1996) found that the time lags decrease with in- 
creasing energy, consistent with our results. But their time lags 
were all negative, while we also find positive timelags. Positive 
time lags were also measured by |Ravasio et al. (2004)| and by 
Brinkmann et al. (2005) in the energy range 0.6-10 keV using 
XMM data. But the statistical significance of our time lags is 
marginal, so one can question their reality. The importance of 
this question and the ample number of lightcurves at different 
energies measured during our multiwavelength campaign would 
deserve a more detailed analysis, which is, however, not the 
scope of this paper. 

The time lags must be produced by energy-dependent mech- 
anisms, for instance, particle cooling, and acceleration. The 
higher-energy particles cool faster and accelerate slower than 
the lower-energy particles, therefore the negative timelag can be 
used to estimate the magnetic-field strength from the measured 
negative time lag. The theory of this process was developed by 
Kazanas et al.[l998 Based on the results of this theory Chiapetti 
et al. d 19991 1 derived the equation (O which provides an estimate 
of the magnetic field B: 



300- (I±*)' 
V -5' 



1 



\ v 



Fig. 9. Intensity correlation between OMC and IBIS data. 
The count-rate range of the ISGRI count rates was split into 
arbitrarily-chosen intervals given by the error bars. 



For determining the frequencies vrj and v\ the mean energy 
for the energy intervals given above were calculated using the 
power laws derived from the fit to the data (see Tables [5] and [6j. 
Inserting those, together with the other parameters from above, 
one gets a magnetic-field strength of (0.08 + 0.03) G. Similar 
values (0.12 G) were found by Kino et al. (2002) using the ob- 
servables of Mrk 421 from Kataoka (2000), by Konopelko et al. 
(120031 (0.1 G) and by Giebels et al. d2007b (-0.1 G). 

But the value of the Doppler factor is highly uncertain. 
Rebillot et al. (2006) investigated two models with 6 = 50 and 
1000 (in their model a higher Doppler factor corresponds to a 
smaller emission region). If we adopt these values, we find val- 
ues of (0.06 + 0.02) G and (0.02 ± 0.007) G, respectively, for the 
magnetic-field strength B. Values of the Doppler factor between 
50 and 100 seem to be suggested by the data (see Rrawczynski 
et al. 120011 and Mastichiadis and Kirk I1997I ). However, such 
large Doppler factors are not observed with VLBI observa- 
tions (e.g. Pin er and Edw ards 2005}, so this might indicate a ve- 



~ c 300 V 

~ ' vi -6> 



(6) 



locity structure in the jet (|Georganopoulos and Kazanas 2003 
Ghisellini]itaL2005J. 



3.5. Correlation analysis 

We investigated whether the X-ray data of ISGRI correlate 
with the optical data of OMC. For this, investigation time in- 
tervals were specified during which the background-subtracted 
ISGRI count rates of Mrk 421 were lying in well-defined, but 
arbitrarily-chosen count-rate intervals. The average ISGRI count 
rates and the average optical fluxes of OMC were calculated 
for these time intervals. The result of this correlation analysis 
is shown in Figure [9] Although the scatter of the results is quite 
large, a slight increase in the optical flux with increasing ISGRI 
count rate is detected. If one fits a straight line to the data, this 
dependence becomes more evident. The functional dependence 
of the OMC flux on the ISGRI rate is 



where 6 is the Doppler factor, z the redshift, and Vo and v\ are 
the frequencies of the corresponding energy intervals (in units of 
10 17 Hz) at which the time lag t (t in s) has been measured. 



= 2.5-10~ 14 + 1.96-10- 



(7) 
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Fig. 10. Result of the fit of a broken power law to the JEM-X 
and ISGRI data in the quiescent state. The top panel shows the 
spectral data with the fit function, the bottom panel the residuals 
of the fit. 
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Fig. 11. Result of the fit of a broken power law to the JEM-X and 
ISGRI data in the active state. The top panel shows the spectral 
data with the fit function, the bottom panel shows the residuals 
of the fit. 



4. Results of spectral analysis 

The results of a spectral analysis of the data of Mrk 421 at X- 
rays are presented in this section. First, we investigated how the 
energy spectrum evolves when going from the quiescent into the 
active state. Secondly, we calculated hardness ratios and looked 
for possible correlations. And third, we combined all data across 
the frequency range and constructed an energy-density spectrum 
v F„. 

4. 1 . Spectral fits 

The data from JEM-X and IBIS were collected for the quiescent 
and active states according to the times given in Table [3] Then 
different spectral models were fitted to the quiescent and active 
fluxes using the functions available in XSPEC11. Since ISGRI 
and JEM-X were excellently cross-calibrated by the two respon- 
sible teams, no cross-calibration factor had to be applied. The 
results of this fitting exercise are given in Tables [5] and [6] In ad- 
dition, a spectrum with a log-parabolic function was fitted to the 
data (see Massaro et al. 2004 ). The results of these fits are shown 
in Table 

It turned out that a broken power law of the form 

FiE) = A ii^vV forE - Eb (8) 

F{E) = A .Ef- a .(-£ v fforE>E b (9) 

(with A the normalization constant and with a and ft the low- and 
high-energy spectral indices, respectively) gave the best fit result 
in both cases. (The joint power law gave an equally good fit, but 
it leads to an unrealistic high value of 6.5 for the high-energy 
power-law index /3 and was therefore discarded.) 

The results of the fits for the quiescent and active states are 
shown in Figures [10] and Q~T] respectively. With a reduced y 1 
value of -1.9, the fits are not excellent, but reasonable. The rea- 
son for the imperfect fit is that the fluxes of JEM-X and ISGRI in 
the overlapping energy intervals differ as is seen from the resid- 
uals. The results of the spectral analysis can be summarised as 
follows: 

- a broken power law fits the data best in both cases (quiescent 
and active) 



Table 8. Results of fits to the Swift-BAT data with a simple 
power law for both states for a reference energy of 1 keV. 



state of source 




normalization constant 




a 




[keV- 1 cm- 2 s" 1 ] 






quiescent 


0.681 


0.53 ± 0.2 


2.46 


±0.11 


active 


1.42 


1.3 ±0.4 


2.49 


±0.14 



- the break energy Eb has a value of ~43 keV 

- the index a is larger than 2 

- the spectral parameters do not change significantly with the 
state (although a slight spectral hardening of the low-energy 
index a is observed) 

The third point indicates that we have only measured the de- 
clining part of the X-ray spectrum and can therefore not deter- 
mine the maximum or minimum of the energy-density spectrum 
v F y . We can only say that the maximum is below the low-energy 
limit of JEM-X of 3 keV. It should be noted here that peak val- 
ues up to 5.5 keV were measured with BeppoSAX (Massaro et 
al. I20041 I. The spectrum also falls off continuously to the upper 
energy limit of 200 keV, so a possible minimum of the v F v spec- 
trum must lie at an energy greater than 200 keV. 

The data from Swift-BAT were also spectrally analysed for 
the quiescent and the active state. They could be best fit with a 
simple power law. The results are summarised in Table[8] Within 
the errors, the spectral indices are identical for the quiescent and 
active states. This is consistent with the results of the data of 
Tables [5] and [6] for the joint power law where the spectral index 
also does not change with intensity. However, the spectral in- 
dex obtained from the Swift-BAT data is with ~2.5 somewhat 
larger than the one obtained from the JEM-X and IBIS/ISGRI 
data (~2.3) but consistent with our result within 3cr. 

The intensity of the Swift-BAT data between the quiescent 
and active states changes by a factor of 2.45 + 1.2. This is a bit 
higher than the values derived from the results of Tables|5]and|6] 
but still compatible within the error. 

4.2. X-Ray luminosity 

By integrating the broken power-law spectrum over the ener- 
gies 3-200 keV and multiplying with the surface of a sphere 
one obtains the isotropic X-ray luminosity in this energy range. 
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Table 5. Model fits to the JEM-X and ISGRI data of the quies- 
cent state. The reference energy is 1 keV. 



Model" 




normalization constant [keV 1 cm 1 s 1 ] 


a b 


PL 


E b [keV£ 


power law (PL) 


2.17 


(0.378 ± 0.004) 


2.318 + 0.006 






joint PL 


1.85 


(9.3 ± 3) • 10- 6 


2.298 ± 0.007 


6.5 ± 1.0 


149 ±11 


broken PL 


1.83 


(0.365 ± 0.005) 


2.300 ± 0.006 


3.0 ± 0.2 


45 ±4 


Band model 


2.2 


(2.8 ± 5) • 10~ 5 


2.0 ± 0.2 


2.32 ± 0.01 


13 ±2 


PL with exponential cutoff 


2.17 


(0.355 ± 0.006) 


2.27 ± 0.01 




311 ± 77 



" The definitions of the mathematical forms of the used functions are 
contained in the XSPEC1 1 manual [with the exception of the joint PL 
which has the form (f-)~°7(l + {§-f~ a )l 

* a is the spectral index at energies below 

c /? is the spectral index at energies above E;, 

d Et is the break energy 



Table 6. Model fits to the JEM-X and ISGRI data of the active state. The reference energy is 1 keV. The parameters and functions 
have similar meaning to the one in Table [5] 



Model 




normalization constant [keV 1 cm 1 s 1 ] 




a 




P 


E b [keV] 


power law (PL) 


3.67 


(0.530 ± 


0.005) 


2.164 


± 


0.004 






joint PL 


1.85 


(9.9 ± 3) 


• io- 6 


2.298 


± 


0.007 


6.65 ± 1.1 


145 ± 10 


broken PL 


1.89 


(0.490 ± 


0.005) 


2.12 


± 


0.05 


2.90 ± 0.08 


41 ±2 


Band model 


3.85 


(9 ±7) 


io- 5 


1.81 


± 


0.09 


2.20 ± 0.02 


23 ± 3 


PL with exponential cutoff 


2.36 


(0.530 ± 


0.005) 


2.04 


± 


0.01 




132.0 ±0.3 



Table 7. Results of fits to the JEM-X and ISGRI data with a log-parabolic function [= K ■ (E/E l y~ b ' log( - E/El ) ] for both states for a 
reference energy of Ei = 40 keV. 



state of 


Xred 


normalization constant 


a 




b 




source 


[keV- 1 cm" 2 s~'] 










quiescent 


2.23 


(7.21 ± 0.2) • IQ-" 


2.35 ± 0.02 


(1.1 


±0.6) 


10"^ 


active 


2.91 


(1.77 ±0.03) • 10~ 4 


2.28 ± 0.01 


(5.0 


±0.5) 


io- 2 



Using the parameters of Tables[5]and[6]one finds luminosities of 
1.75-10 45 erg/s and 3.73-10 45 erg/s for the quiet and active states, 
respectively. This is much lower than the Eddington limit, which 
is in the range of (2.4 - 9.6)- 10 46 erg/s. But it should be noted 
that the emission is actually beamed and that the luminosity is 
thus lower than given above. 

4.3. Hardness ratios 

We calculated the hardness ratio (HR) of the two IBIS/ISGRI en- 
ergy bands 20 - 40 keV (L-band) and 40- 100 keV (H-band) HR 
= H-band rate / L-band rate. It is plotted in Figure[T2]as a func- 
tion of the time and in Figures [L3l and [T4l as functions of the L- 
band and the H-band intensities. The hardness ratio in Figure[T2l 
is fairly constant in time and does not follow the lightcurves of 
Figures [5] and [6] so we do not observe an evolution in the hard- 
ness ratio with intensity. When inspecting Figures [l3]-[T4j one 
notes a remarkable difference: whereas no correlation is obvious 
when one plots the ratios as a function of the L-band count rate 
(Figure [13) , a correlation for intensities up to ~2 counts/s seems 
to be present when the plot is done as a function of the H-band 
count rate (Figure [141 . But this might be because there are no 
events with <1 counts/s in Figure [L3l 

In Figure[l5]the hardness ratio is shown for two energy bands 
[(l-10keV)/(0.3-l keV)] of the Swift-XRT. A comparison with 
the lightcurve of Figure [3] shows that the hardness ratio follows 
the lightcurve closely. This means that we have observed an 
intensity-hardness-ratio correlation; i. e., the brighter the source 
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Fig. 12. Hardness ratio H/L as a function of time for the 
IBIS/ISGRI data. 



the harder the spectrum. But this is obviously not valid at all 
times. The last peak in Figure [3] (at MJD m 53910.5) is only 
marginally visible in Figure [15] One would expect a hardness 
ratio of ~2 (similar to the hardness ratio observed at MJD » 
53909.5), but actually it is only about 1.5. This means that the 
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Fig. 13. Hardness ratio H/L as a function of the L-band (20 
40 keV) intensity for the IBIS/ISGRI data. 
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Fig. 16. Hardness ratio of the Swift-XRT data as a function of 
the XRT intensity. 
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Fig. 14. Hardness ratio H/L as a function of the H-band (40 - 
100 keV) intensity for the IBIS/ISGRI data. 
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Fig. 15. Hardness ratio of the Swift-XRT data as a function of 
the time. 



emission process must have changed significantly on a time scale 
of <1 day. We obviously observed a clear hard-to-soft evolution. 

The overall hardness-ratio correlation is also visible when 
one plots the Swift-XRT hardness ratio as a function of the 
Swift-XRT count rate as in Figure[l6] Although the positive cor- 
relation with the rate is not striking it is clearly recognisable. The 
large scatter around this correlation is probably due to a short- 
time scale and overlapping variability. 

4.4. Multiwavelength spectrum v F v 

The data from all observations (with the exception of the JEM-X 
and ISGRI data) were averaged over the observation time span, 
corrected for extinction at optical wavelengths using the extinc- 
tion value Eb-v — 0.03 of Burstein and Heiles (1982) and the for- 
malism of Seaton (1979), converted to the same unit [erg/(cm 2 
s)], and then plotted in an energy-density spectrum (v F v spec- 
trum). This multiwavelength spectrum is shown in Figure [T71 
Unfortunately the number of data points is sparse, especially at 
energies above -500 keV where only two data points (measured 
at different times) at TeV energies exist. It should also be noted 
that the data used to create the average high state and low state 
SED data points shown in Figure [17] were extracted from the 
same time windows for all instruments; however, the data are 
not strictly simultaneous. Due to this non-simultaneous data, it 
is possible that small time-scale flaring could cause a system- 
atic offset of the SED from one energy band to another. This 
could explain the discrepancy between the Swift-XRT data and 
the JEM-X and IBIS data. 

The data are compared with the theoretical models of 
IMaraschi efaLl (1999) , Kino et al. (120021 ) and Blazejowski et al. 
(2005 ). The models of Maras chl et al.l (1999) and Blazejowski et 
al. (120051) were adjusted to the flux measured by IBIS in the qui- 
escent state by applying a factor of 3.2. Both models predict the 
actual flux measured by Whipple at TeV energies quite well, but 
seem to fail to fit the optical data and the radio flux. On the other 
hand, the model of Kino et al. (2002) fits the radio data quite 
well, but has problems fitting the X-ray and TeV data, together 
with the radio data. It is remarkable that the JEM-X spectrum is 
flatter than predicted by the models, even though their shape dif- 
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Fig. 17. The measured energy-density spectrum of Mrk 421 is compared with the three theoretical models of Maraschi et al. (1999), 
Kino et al. (120021 1, and Blazejowski et al. (2005 ). Please note that the data and the three models represent the source in different 
states and therefore reflect its spectral variability. 



fers significantly. The three models also do not predict the break 
at an energy of 43 keV. 

The best agreement between data and models is found for the 
model of Blazejowski et al. (2005 ). However, it should be noted 
here that the models were fitted to data measured at different 
times, so it is not surprising that they do not fit our data well. 
In contrast, they show that not only the intensity, but also the 
spectral shape of Mrk 421 is extremely variable. 



4.5. Model fits 

We also tried to fit a theoretical emission model to the data of 
the energy-density spectrum of Figure [17] The model curves 
fitting the data shown in Figure [18] best were obtained using a 
one-zone synchrotron self-Compton model (SSC) including the 
full Klein-Nishina cross section for inverse Compton scattering 
( IJones et a l. 1968 Blumenth al et al. 19701 1 . The model assumes 
a spherical blob with constant injection of non-thermal electrons 
and a constant escape rate for electrons and photons. The model 
details are described in Ruger et al.|(2007). 



In this calculation the electron distribution is assumed to be 
constanst (stationary) and is represented by a smoothed broken 



power law combined with an exponential cut-off as given in 
|Tavecchio et al. (200 1)| 



N(y) = Ky- 



Jb 



exp 



(10) 



with the differential electron number density N, Lorentz factor y, 
break energy £), = yt>mc 2 , spectral indices s\ and S2, normaliza- 
tion factor K, and cut-off energy y max . This phenomenological 
Ansatz for the spectrum is justified by self-consistent simula- 
tions (see Ruger et al. 2007). 

With this electron distribution, the radiative transfer equa- 
tion was solved numerically in the comoving frame and the ob- 
served spectrum obtained by applying the beaming effect using 
the Doppler-factor 6. Specifically, the transformation for the in- 
tensity is 



*v,beamed 



6*1, 



y, comoving frame- 



(ID 



In addition, one can show that, for dominating synchrotron 
losses, the second spectral index can be rewritten as si — s\ + 1, 
so one less parameter has to be determined. A first estimation 
of the parameter set is obtained by using the equations given 
in Kataoka et al. (2000) Small variations lead to the best-fit pa- 
rameter set presented in Table [9] 
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Fig. 18. The result of the model fits to the measured energy-density spectrum of Mrk 421. 



Table 9. Best-fit parameters for high-state and low-state inten- 
sity when a one-zone synchrotron self-Compton model is fit to 
the spectrum. 



parameter 


high-state 


low-state 


Jb 


4- 10 s 


3 • 10 5 


Ymax 


1.2- 10 6 


7 ■ 10 5 


S\ 


2.2 


2.2 


Si 


3.2 


3.2 


B[G] 


0.1 


0.1 


K [cm" 3 ] 


1.6- 10 4 


1.6- 10 4 


6 


15 


15 



From this table one can see that both parameter sets differ 
only in the break energy j/, and the cut-off energy y max . The latter 
energy represents the efficiency of the acceleration mechanism. 
Since yt, results from the emerging balance between cooling and 
acceleration of the electrons, a more efficient acceleration mech- 
anism leads to both a higher y max and a higher yt,. Therefore we 
conclude that the variability of the emission can be due to vary- 
ing efficiency of the acceleration. 

It should be noted, however, that a change in the spec- 
trum of the injected electrons without changing the environ- 
mental parameters like the magnetic fields and/or the Doppler 
factor is difficult to understand. But the acceleration of par- 
ticles can be influenced without changing the environmen- 



tal parameters by taking turbulence into account which can 
have a dramatic influence on the particle acceleration; on the 
one hand, enhanced turbulence can change the spectral in- 
dices dVainio & Schlickeiser 19981 and Vainio & Spanier 2005 ), 
and on the other, the decay length of turbulence can 
change the maximum energy dVainio & Schlickeise r 19981 and 



Vainio & Spanier 2005 1. However these effects are not covered 
in our model, just like the geometry of the shock and the shock 
thickness. 



The fit presented here is not unique (see discussion in 
ITavecchio et al. 19981 1 and the observed change in spectral en- 
ergy distribution can also be modelled by slightly varying 
the magnetic field strength and Doppler factor. The changes 
in these parameters would lead to change in electron spec- 
trum, and such detailed modelling is beyond the scope of this 
paper. The Doppler factor and magnetic-field strength mod- 
elled here are similar to those derived for other TeV-emitting 
blazars. However, very long baseline interferometric obser- 
vations show that such high Doppler factors would require 
very extreme assumptions for the viewing angle (0.015 de- 
grees) and the Lorenz factor (25) (IPiner a nd Edwards 2005). 
Therefore adopting models where the jet has a velocity struc- 
ture (Georganopoulos and Ka zanas 2003|IGhisellini et al. 2 005) 
might be more feasible, but such modelling is beyond the scope 
of the paper. 
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5. Discussion and conclusions 

The blazar Mrk 421 underwent an active phase in May/June 
2006. INTEGRAL was therefore reoriented to observe this 
source for 12 days with the instruments OMC, JEM-X, and IBIS. 
Simultaneously the source was also observed by the Metsahovi 
and VLBA radiotelescopes, by the KVA telescope in the optical 
waveband, by RXTE and SWIFT at X-rays, and by Whipple at 
TeV energies. In X-rays, several strong flares were observed not 
seen at lower energies. Unfortunately the data at TeV energies 
are too sparse to allow such a conclusion for the high energies. In 
this context it is especially striking that the flares are not visible 
in the optical lightcurves, since it is thought that the X-rays and 
optical photons are produced by the same electron population 
via synchrotron emission. One would therefore expect similar 
lightcurves; but when fitting a theoretical one-zone synchrotron 
self-Compton model to the data, it is found that flares seen in 
X-rays may not necessarily be visible at lower frequencies. 

However, other possibilities also exist for explaining the ob- 
served behaviour, rather than varying the efficiency of the accel- 
eration. For example, it could be that the optical photons origi- 
nate in larger emission regions and that the flares are thus occur- 
ring on longer time scales and are therefore smeared out. It could 
also be that the intensity fluctuations seen in X-ray are damped 
in the optical band. This view could be supported by the ob- 
served intensity correlation between OMC and IBIS in Figure[9] 
Whereas the IBIS rate varies by a factor of ~5, the variation in 
the optical flux is only a few percent, so the high variations seen 
in X-rays are strongly damped in the optical range, hence barely 
visible. This could be achieved, for example, if the optical pho- 
tons were obscured by some intervening clouds of matter that 
are transparent for the X-rays but not for the optical photons. 

The spectral analysis of the JEM-X and IBIS/ISGRI data 
gave a surprising result: the spectral index a of the power law 
for the two states does not change much (from 2. 12 to 2.3) at en- 
ergies below 40 keV, although the intensity changes by factors 
up to 5. In the energy range 2-10 keV, Takah ashi et al. (1 996) 
found a clockwise evolution of the spectral index as a function 
of the intensity. But in their case the range of the change of the 
spectral index is nearly the same as reported by us (from 2.3 
to 2.52), although the intensity only changes by a factor of 2. 
This shows that the observed slight hardening occurs predomi- 
nantly at low intensities. This "saturation effect" is clearly seen 
in Figure [14] 

Another remarkable result is that the hardness ratio of 
the Swift-XRT data is strongly varying with intensity (see 
Figure [TBI , whereas the hardness ratio of the IBIS/ISGRI data 
is more or less constant (see Figure [L2l. This again shows that 
there is a conspicuous difference in the emission processes of 
the low- and high-energy X-rays. At energies <10 keV, the spec- 
trum becomes harder with increasing intensity. This is not the 
case for X-rays with energies >20 keV. In the standard model 
in which the X-rays are produced by synchrotron radiation from 
the same electron population, this different behaviour is difficult 
to explain. 

The picture becomes even more complex when one com- 
pares the lightcurves of Figure [3] to Figure [6] The two flares at 
the end of the observation interval (around MJD 53910) show a 
noticeable spectral difference: the first peak is much harder than 
the second one! If one assumes that each flare is produced by 
a relativistically-moving blob of electrons, the physical charac- 
teristics of each electron population and their environment (i. e. 
magnetic fields etc.) must be different, so the picture of jet emis- 
sion is very complex and the extraction of the relevant informa- 



tion from the data requires a detailed and profound modelling of 
the emission processes. 
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